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It is investigated how a spatial quenched disorder modifies the dynamics of cou- 
pled map lattices. The disorder is introduced via the presence or absence of coupling 
terms among lattice sites. Two nonlinear maps have been considered embodying 
two paradigmatic dynamics. The Miller and Huse map can be associated with an 
Ising-like dynamics, whereas the logistic coupled maps is a prototype of a non trivial 
collective dynamics. Various indicators quantifying the overall behavior, demon- 
strates that even a small amount of spatial disorder is capable to alter the dynamics 
found for purely ordered cases. 
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I. INTRODUCTION 

Far-from-equilibrium systems, in which the simultaneous presence of non-hnear terms and 
spatial extent is an essential feature, involve many important physical or chemical phe- 
nomena as varied as Raylei^h-Benard convection in large horizontal boxes chemical 
react ion- diffusion systems 2'], colonies of micro-organisms and fibrillating heart tissue 
0]. In all these examples, the spatio-temporal chaos arising from the coupling of a large 
number of degree of freedoms, poses challenging theoretical questions which has so far eluded 
a universally accepted description. Progresses have often been limited to models contain- 
ing the essential ingredients of typical extended systems at the crudest possible level. In 
this respect, "simple" partial differential equations such as, the Kuramoto-Shivashinsky, 
Ginzburg-Landau and Swift-Hohenberg equations, have proven to be useful models to iden- 
tify main characteristics such as symmetry breaking, phase transition, phase dynamics, 
spatial chaos or defects. From a numerical point of view, however, these approaches, albeit 
successful, are rather demanding in terms of computational time. Coupled map lattices 
(CML) p, where time and space are discretized from the outset, provide an alternative 
computational fast tool to study dynamical processes in spatially distributed systems. As 
the coupling with nearest-neighbours mimics the diffusion process, CML may be reckoned as 
a non-linear version of the diffusion process on a lattice. As a matter of fact, it is the subtle 
balance between diffusion and nonlinear effects that leads to the patterns and non-trivial 
time behaviors found in CML. 

Pure diffusion on lattices has a long and venerable tradition in statistical mechanics P]. It 
was originally motivated by solid-state applications as toy models for particle transport in 
ordered and disordered solids, but it has soon grown into a field in its own right. Interesting 
issues arise when a quenched disorder is introduced to model the fact that real life diffusion 
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for a review). The 



occurs in disordered media such as a porous medium (see e.g. 
presence of disorder drastically modifies the diffusion properties, sometime in a somewhat 
unexpected way Depending on its type and its strength, the disorder either inhibits 
diffusion (subdiffusion process) or enhances it by creating priviledged paths (superdiffusion). 
Motivated by this scenario, the natural question arises of what would happen for CML 
in the presence of a quenched disorder, where three different mechanisms (non-linearities, 
diffusion and disorder) compete against one-another. This issue has already been addressed 



for various non-linear systems |8|. Stochastically induced synchronization driven by the 
presence of additive noise has also been observed in spatially extended systems 

Clearly the presence of disorder induces a translational invariance breaking which gen- 
erates a particular case of a system with a large number of degrees of freedom which are 
non equivalent. In analogy to what is observed in other problems of statistical physics, new 
phenomena may arise. As a matter of fact, aging, dynamical phase transition, anomalous 
transport have been observed in one-dimensional maps 

In the present paper, we focus on two-dimensional CML. Two such systems are investigated. 
On the regular lattice, the first one displays phase transitions ("ferromagnetic ordering") 
and, the second one, a non-trivial global dynamics Disorder is introduced at the simplest 
possible level, namely a prescribed fraction of bonds are randomly disconnected. This affects 
the phase diagram of the system : even a tiny fraction of disconnected bonds markedly 
modifies or destroys altogether, the dynamics found in the ordered lattice case. 
The paper is organized as follows. In Sect. lUl we briefly review some results on diffusion on 
a lattice and nonlinear coupled maps, while Sect lllll introduces the way disorder is imple- 
mented in the lattice. Sect JIVI details the numerical results for both CML models studied 
and Sect|V]is devoted to some concluding remarks. 

II. DIFFUSION AND COUPLED MAPS ON A LATTICE 

The most general master equation governing the linear diffusion of particles on a lattice 
reads 

dtP.At) = Yl l^-yP-oyit) - Wy,P,,,{t)] (1) 

where P^^r^lt) stands for the probability of being at site x G Z'^ of a hypercubic lattice of 
dimension d, having started at site Xq at the initial time to = 0. Wxy denotes the transition 
rate for a jump y ^ x and is different from zero only for nearest-neighbours y of site x (the 
notation y{x) in the sum reflects this constraint). A time discretized version (in time units 
such that At = 1) reads 

PxAt+'^) = (1 -ex)i'a=„x(t) + Xl^-J/^-o?/W where = ^ W^j,, (2) 

y{x) y{x) 
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If only positive symmetric transition rates < W^y = Wyx are allowed and ^^^(a,) W^y = 1 
(for all x), then one easily verifies that (a) < Px^xit) < 1 and (b) the total probability 
is conserved. The physical process associated to such a mathematical framework can be 
described as follows : at each time t a particle at x can either move to a nearest-neighbour 
y{x) with rate Wyx or stay at x with probability I — ex- Quantities Wxy are generally random 
with the proviso that condition < e^^ < 1 should hold true for any x. Standard diffusion 
corresponds to Wxy = 1/z where z = 2d denotes the coordination number of a d-dimensional 
hypercubic lattice. Introducing a quenched disorder in the transition rates Wxy changes the 
classical diffusion behavior into a sub-diffusive or super- diffusive one jrll. Such a toy model 
may be used to understand a variety of complex systems such as optical excitation, polymers 
and diffusion-limited binary reaction ^. 

Other interesting dynamical features arise when non-linear terms are introduced in the local 
dynamics. Let Cx{t) be a real variable defined at point x and time t. This quantity belongs 
to the interval [—1, 1] and satisfies 

Cx{t + 1) = (1 - zg)F [Cxit)] +gY,F [^J/ W] (3) 

in which F{C) stands for a non-linear homorphism in [—1, 1]. Quantity Cx{t) may represent 
the concentration of a binary mixture, the values —1 and -|-1 corresponding to a pure phase of 
the first or second component respectively. The above equation resembles equation Q with 
Wxy = g for all nearest-neighbours but with a non-linear map F{C). Coupling g between 
neighbour sites is constrained to satisfy < g < l/2(i so that, if the initial concentration is 
in the range [—1, 1], it remains within the same interval for all times. The first r.h.s factor 
in equation (jH)) accounts for the "self-interaction" term, that is the coupling of each site 
with itself. As g 0~^, each site evolves independently of the other ("high temperature 
regime"). At g = l/{2d + 1), each site is coupled with all its nearest-neighbours as well as 
with itself in an identical manner ("democratic coupling"). Finally as g ^ gmax = l/2(i, 
the self-interaction term vanishes and all nearest-neighbour sites are strongly and equally 
coupled ("low temperature regime"). For a two-dimensional lattice, the democratic coupling 
and the low-temperature regime correspond to g = 0.2 and g = 0.25 respectively. 
In the present work, we consider a finite size lattice consisting of L sites along each spatial 
direction with standard periodic boundary conditions. Two maps F{C) are studied, each one 



representing a particular class of dynamics. The first one is the Miller and Huse map 
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-2/i/3 -fiC for - 1 < C < -1/3 

fiC for - 1/3 < C < 1/3 (4) 

2/i/3 - /iC for 1/3 < C < 1 



where C G [—1, 1]. The map (@]) with /i = 3 was introduced by Miller and Huse 
to mimic a statistical system having two equilibrium states i.e. a Ising-like dynamics. 
This model is appropriate to describe a spinodal decomposition mechanism if variable 
ax{t) = sign [Ca;(t)] = ±1 plays the role of an Ising spin. In that case, a magnetization 
can be defined as follows: 



where = L'^ is the total number of sites. 
The second study concerns the logistic map 



F^{C) = l-ixC^ 0</x<2 ^€[-1,1]. (6) 



This CML possesses a number of interesting dynamical features In particular, the 

spatially averaged quantity 

W) = lpJ2^.{t) (7) 

X 

possesses a macroscopic time-periodic evolution while spatial patterns similar to those found 
for Miller-Huse map are not visible. Moreover the graph of the asymptotic time evolution of 



C(t), with respect to parameter n displays a bifurcation diagram similar to the one observed 
for a unique Logistic map Ql- 

III. COUPLED MAP LATTICE AND DISORDER. 

In close analogy with linear diffusion on a lattice, the effect of a quenched disorder is intro- 
duced on the two CML discussed above. The generalization of Eq.(jni) then reads 

a(t+l) = (l-ejF[a(t)]+ J]W,,F[C,(t)] (8) 
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Paralleling the diffusion case, couplings are partially random and verify for arbitrary x 

Wy, = W^y, < = J2Wyx < I. (9) 

As a consequence, quantities C^(t) always remain within the original interval [—1,1] since 
\C,{t + 1) I < (1 - e,.) + J2 < (1 - + J2^y- = 1- (10) 

y{x) y{x) 

This model includes all the previous cases: 

a) Linear diffusion with disorder (j2)) when F{x) = x. 

b) Nonlinear coupled maps © when W^y = g (and hence = zg). 

In the following, a quenched binary disorder is used satisfying the distribution 

p{w) = (1 — p) S{w — g) + pS{w) (11) 

where p stands for the percentage of connections which are cut. The value p = corresponds 
to the regular case. It should be emphasized that, as bonds are removed, the self-interaction 
term increases (see Eq.(jH))). Disorder may set a competitive mechanism disfavouring the 
balance involving couplings between lattice points, which tends to "synchronize" the maps 
spatially or temporally and the non-linearity of the maps, which tend to destabilize them. 
Disorder is thus expected to act as a further destabilization preventing the onset of an 
"ordered" phase. 

IV. NUMERICAL RESULTS 

Numerical results are obtained on a square lattice {d = 2). The Miller and Huse case is first 
presented in Sect. A,B,C as a prototype of a spinodal decomposition mechanism. Sect. D 
and E are devoted to the logistic map where a global time evolution is observed. 

A. The ordered Miller and Huse phase diagram for /i = 3 and /i = 1.9 

FigG] (top) depicts the Miller- Huse phase diagram for /i = 3 |l^ . If denotes the value of 
the mean-square magnetization, this quantity is raised to two suitable exponents (m^)^ and 
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{L'^m?') ^ , and plotted as a function of the coupling strength g. In their original paper, Miller 
and Huse used the two dimensional Ising critical exponents (3ising = 1/8 and = 7/4 

jisl ]. A successive and more refined work jl^ identified the characteristics exponents values 
/3 ~ 0.111, 7 ~ 1.55 which are close but not identical to those of the two dimensional Ising 
model. Figure ^ was obtained for lattice size L = 128 with a magnetization which was 
computed after T = 2 x 10^ iterations which appear to be large enough to reach a steady 
state, and averaged over 10 — 50 different initial conditions (50 refers to the left branch of the 
phase diagram g < g^, where fluctuations are much higher). As one approaches the critical 
value gc = 0.20(5), a phase transition occurs. This result is in agreement with those of 
which pinned down the transition value at gc = 0.2051(5) on the basis of more extensive 
simulations. FigUl (bottom) provides the same computations for /i = 1.9 and indicates a 
critical value gc = 0.16(8) not previously reported. Here and below the parenthesis indicates 
the degree of uncertainty of the result. 

The fact that the total magnetization is not conserved, suggests that this dynamics could be 
akin to a dynamical Ising model with non-conserved magnetization (model A) jiol. [l7 1 . How- 
ever, the critical coupling, as well as the time evolution (see below) gc changes for different 
fi. This feature is apparent in Fig. |21 where, starting from a random and uncorrelated initial 
condition for g > gdfJ'), two snapshots are depicted for = 3 (top) and fi = 1.9 (bottom) 
at stages t = 100 (left) and t = 1000 (right). 



B. The persistence probability in the case of a regular lattice 

Another interesting quantity, is the persistence probability Ps{t), which measures the fraction 
of spins which have not changed sign up to time t. This concept was first introduced within 
this context in It was shown that, for fi = 1.9, a transition occurs at gc ~ 0.169 

appearing as an abrupt modification in the slope of probability Ps(t) (Fig El). For gc < g, 
this probability takes the form 

Ps{t) ~ t-' with e{g)r^ig-gcr (12) 

Values w ~ 0.20 and gc = 0.16(9) are reported in [l]]]. This latter number is strikingly close 
to the critical value gc = 0.16(8) found in the above section. The transition of the persistence 
probability thus appears to be a consequence of the phase transition of magnetization (and 



not computed in rather than a different transition. 
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C. Miller-Huse coupled map lattice with binary quenched disorder 

Let us now introduce a shght disorder on the Miller and Huse coupled maps. Two small 
percentages of disconnected bonds have been used : p = 0.05 and p = 0.065. Results for 
/i = 3 are reported in Fig. HI Note that exponents of the undiluted case are employed on 
the basis that critical exponents of the Ising model with random binary coupling are found 
to be similar to those computed for a regular lattice Dilution has the effect of shifting 
the critical coupling to higher values. This shift can be qualitatively understood by a 
simple argument. Let denote z{p) = z{l—p) the average effective coordination number with 
dilution p. Within a mean field description, the average coupling strength increases in the 
presence of dilution by a factor 

gcip = 0) I-P ^ ^ 

thus leading to a linear dependence for gdp) for small values of p. Our calculations reproduce 
this linear behavior (Fig. Ej), although with a different numerical value of slope. 
As the maximum possible g is equal to a constant (7max = l/2d = 0.25, the right branch of 
the phase diagram shrinks as p increases, and there will be a critical threshold Pc, such that 
g{Pc) = fl'max, above which a phase transition is no longer observed. Results of FiglHlthen 
suggest Pc ~ 0.11, although a more extensive analysis would be needed for a precise estimate. 
We checked that, for values p > pc, the ferromagnetic transition is indeed not observed. A 
small percentage of disconnected bonds therefore disfavours a collective behavior and can 
even inhibit the transition. A similar trend is mirrored by the persistence probability (Fig. 
in)), which is computed with a dilution p = 0.05 and under the same conditions as in Fig. 01 
This also supports the equivalence between persistence and ferromagnetic transitions. 
Patterns for both /i = 3 and fi = 1.9 above the critical coupling {g > gdf^)) are shown at 
some intermediate time {t = 1000) in Fig. [3 for dilution p = 0.05. These are the analogous 
of those previously obtained under the same conditions in the regular lattice (Fig. E}. If 
patterns are pinned by defects, we expect similar patterns upon a rescaling conserving the 
quantity pL'^. Indeed, the number of "defects" (sites with disconnected bonds) should be 
equal to the fraction p of disconnected bonds times the total surface N = L"^. Figure |H1 
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indicates that this is precisely the case as the patterns obtained at a relatively late stage of 
the evolution (t = 10^) L = 128, p = 0.05, and L = 256, p = 0.0125, L = 512, p = 0.003125 
are quite similar. 



D. The Logistic CML and non-trivial collective behavior 

In the Miller- Huse CML, the emphasis was on collective spatial behavior in the sense of 
ordinary critical phenomena. By contrast, the logistic CML (jH)) is an example of a collective 
time evolution [ll|. Depending on parameter /i G [0, 2], a large variety of complex behaviors 
are observed. Building upon previous analyses on this system with no disorder jll|, an 
extensive analysis is performed both with and without disorder. Our results are in perfect 
agreement with published results [ill !3- addition, we address issues such as transient 
behavior, coupling dependence and competition with disorder, that were not previously 
covered. 



System ^ is first iterated for extremely long time and spatial average C{t) is thereafter 
computed when a steady state is expected to be attained. Note that some numerical tests 
require times up to T = 5 x 10^ which would be prohibitively long for size L = 2048. Hence 
our investigations have been performed with smaller sizes {L = 1024, L = 128 and L = 64). 



This steady state condition is checked for by considering C (t) at several successive large times 
and showing that the same sequence is periodically found. This sequence is then analyzed 
as a function of parameter fi for different couplings g. For democratic coupling g = 0.2, 
a highly regular time evolution arises yielding a periodic doubling bifurcation cascade for 



C{t). Fig. El reports results for two computations: L = 128 with T = 10^, and L = 1024 
with T = 10^. Both results nicely superimpose and are compatible with those of 
up to fluctuations of order 1/VL (in these work, a size L = 2048 was used). This suggests 
that the CML behaves in the thermodynamic limit and that steady state has been reached 
in both cases. For other couplings, the observed behavior can be more involved. Indeed 
identical computations were repeated for weaker couplings {g = 0.1 in Fig. E)). For L = 128, 
the same final time is used, and quite a few "singular" points (i.e. outside the bifurcation 
diagram) are now distinguishable. For size L = 1024, computations have been performed 
up to times as large as 10^ and the system is always found in a transient state (see last 
picture of FiglHl), in contrast with the L = 128 case. As in the uncoupled case g = 0, values 
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of Cx{t) are almost uniformly distributed with the exception of some particular windows 
characteristic of the single map, it is then plausible to expect longer transient regimes and a 
stronger dependence on the initial conditions in the weak coupling regime [g < 0.2). Indeed 
Fig. El shows snapshots of spins ax{t) (left) and fields Cx{t) (right) at two extremely long 
times (T and T + 1 with T = 10^) for identical parameters {L = 64, g = 0.1, fi = 1.66). The 
two situations only differs by their random initial conditions. For the first initial conditions 
(first two rows), an oscillation is observed between two different states corresponding to the 
two branches of Figure IHl This is confirmed by the corresponding probability distribution 
function (PDF) in which a well defined bimodal character can be identified(Fig. ^2 (a)). 
The second random initial condition leads to a different pattern (see the last two rows of 
Fig llOj) . The lattice is now divided into two regions. In each region, the system is behaving 
as before but there exists a relative phase lag between the two " sub-lattices" . This feature 
is also visible in the PDF, of Fig^J (b) where the new bimodal PDF is a linear combination 



of the previous two distributions. As a consequence, the value of C{t) is shifted towards 
"singular" points reported in FiglHl These "frozen" states have extremely long and possibly 
infinite relaxation times. 

Upon identifying a suitable order parameter and monitoring the variation of such a parame- 
ter with coupling g, one finds a transition from a "high temperature" to a "low temperature" 
phase, akin to that found for the Miller-Huse map. A proper definition of the order param- 
eter is, however, a delicate matter. Depending on the value of the coupling g, the spatial 



average C{t) is a time constant or it follows a periodic sequence. Hence, by computing 



Cit+1)-C{t) 



C{t+1) + C{t) 



with Ca,e{t) = ^ (14) 



Cave (^) 

for sufficiently long times, it is expected to get value of ip{t oo) close to for small g, 
and a non-zero value for g > gc (iSf. An abrupt change is indeed observed (Fig. IT^ for a 
given gdfJ')- For instance, one obtains (^c ~ 0.1 for /i = 1.7. More extensive and systematic 
calculations on this point are left for a future study. 



E. Effect of a binary quenched disorder on the logistic coupled map lattice. 

Let us now consider the coupled logistic maps in the presence of binary quenched disorder. 
We find that even a small amount of disorder is sufficient to modify the "conventional" 
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pattern reported on the regular lattice, and this is true irrespective of sizes L. Moreover, 
long transients are now present even in the democratic coupling g = 0.2. For instance we 
have explicitly checked that, for L = 1024, T = 10^ is no longer sufficient to reach the steady 
state. We hence consider size L = 128 for which a sufficiently long time (T = 4.5 x 10^) 
can be attained to ensure steady state conditions. Spatial average C{t) is shown in figure UHl 
versus parameter /i, for coupling g = 0.2 and two weak dilutions p = 0.05 and p = 0.1. 
Once a small fraction of the total number of bonds is missing, the formation of the period 
doubling sequence is shifted to the left : disorder thus tends to inhibit the formation of 
the bifurcation structure. An additional increase in the dilution p gives rise to a further 
shift to the left, and eventually to destruction of the collective dynamics for a sufficiently 
high p = pc- This conclusion is supported by further calculations (not shown) on a smaller 
size (L = 64) where a systematic trend to the left is observed as p increases. However, the 
large fluctuations present for such a small size prevent us from a quantitative measure of 
the critical threshold. 



F. The Lyapunov spectrum 

I — I 

As suggested in ^20] , the Lyapunov spectrum is a priori another useful quantity to analyze 
the collective behavior of CLM. The computation is performed by recasting Eq.(jH} in the 
following form: 

a(t + l) = ^M,j,F[C,(t)] (15) 
y 

where the matrix M^y stands for 

M.y {{W}) = (1 - e^)6^y + W,y5\,.y\^l (16) 

Upon considering a variation 5C{t) of a given solution C°(t), one obtains 

5a(t + 1) = E ^^y^' [C'J(t)] 5Cy{t) (17) 

y 

The map enters into equation (fTTj) only through the value The Lyapunov exponents 

are computed by taking into account the time evolution of an initially small variation 6Cx{t = 
0). This procedure is standard and outlined in j2l|. To compute the lower part of the 
spectrum with sufficient precision, the time window for the evolution is kept sufficiently 
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small and the number of time intervals for normalization is increased accordingly. This 
procedure has been applied for both the Miller-Huse and the logistic map. A cross-checking 
for the good performance of the algorithm can be carried out by computing the exact 
eigenvalues for the uncoupled case g = and the diffusion linear map . In this latter case, 
the eigenvalues of matrix M^y are directly connected to the Lyapunov exponents. 
In the regular case, the Lyapunov exponent spectrum has already been carried out for 
the Miller-Huse map at yU = 3 [2^. In that work, it was suggested that the Lyapunov 
spectrum is not sensitive to the presence of the "ferromagnetic ordering" discussed earlier. 
The independence of the spectrum on the critical coupling gc was also related to the fact 
that chaotic fluctuations, as measured by the Lyapunov spectrum, have short range and 
thus decouple from the the long-range spatial effect giving rise to the ferromagnetic ordering 
at the critical coupling g^. This result is confirmed by our computations which were carried 
out at various sizes (from L = 24 as in ^2^ to L = 256). Figure IT^ plots the first 100 
most unstable exponents (this represents a small fraction of the spectrum whose size is 
V = L"^ = 16384 for L = 128) versus the normalized index x = z/L^ for two values of the 
coupling (one below and the other above the critical coupling gc ~ 0.205). As in [20], no sign 
of transition can be traced in the results. The spectrum is continuous starting with a value 
which is well below the unique exponent of the uncoupled case ~ 0.9. The independence 
of these results for the upper spectrum on the size L of the system has been verified. The 
effect of the dilution on this spectrum has been analyzed in the presence of different degrees 
p of dilution. Fig. HH displays the results for p = 0.1 showing that disorder lowers the values 
of the Lyapunov spectrum. 

For the logistic map, we are not aware of any previous calculations even in the regular case. 
Fig. El depicts our results both in the undiluted and the diluted case in the case of a 
democratic coupling g = 0.2 and for two different values fi = 1.85 and fi = 1.65. For the 
first /i the spatial average C{t) is a time constant, whereas for the second value of /x, C{t) 
is time periodic. Nevertheless, no qualitavive difference is found in the Lyapunov spectrum 
which is simply shifted to higher values as /i increases. Unlike in the Miller-Huse case, the 
presence of disorder appears in the Lyapunov spectrum in the form of a gap separating the 
first few Lyapunov exponents from the rest of the spectrum (Fig. [T3|) . Again, the effect of 
changing the value of the parameter is a shifting of the numerical value of the exponents, 
in agreement with the undiluted case. On recalling the analogy between the Lyapunov 
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spectrum and the problem of the Schodinger equation on a lattice [22], the emergence of 
the gap in the lowest part of the Lyapunov spectrum in the presence of "defects", can be 
regarded as the analog of a "localization" mechanism. On the other hand this analogy is 
certainly limited, or it is at least incomplete, since it does not apply to the Miller and Huse 
case discussed earlier. 



V. CONCLUSION 



We have shown that the phenomenology of spatio-temporal nonlinear systems may be signifi- 
cantly changed by a small amount of quenched disorder. This feature has been demonstrated 
on two-dimensional CML with periodic boundary conditions. Two paradigmatic cases have 
been investigated. The first one is the Miller and Huse case, which is known to display a 
phase transition in the strong coupling (low temperature) phase. The transition is found 
to be delayed and even inhibited by a small amount of binary noise. A percentage of dis- 
connected bonds as small as 10% is sufficient to destroy the ordered phase. This feature is 
mirrored by an analogous transition on the persistence probability. The second paradigmatic 
case concerns the logistic map which exhibits a remarkable time evolution of the spatial av- 
erage C{t) on a regular lattice. In this case, issues regarding the possibility of long transient 
behaviors have been addressed. As it turns out, in the non-democratic case, the system 
might get trapped into a metastable state away from the expected bifurcation branch. In 
addition, the phase diagram C(t) versus /i is markedly affected (even in the democratic case) 
by a small percentage of dilution (e.g. 5%). 

This work is believed to unravel some aspects of the interplay between quenched disorder 
and non-linear processes for cases which were not previously studied. In particular, it 
indicates that the lattice structure affects in a major way the system dynamics. Along 
this line of thought, it would be worth modifying the nature of the disorder to understand 
the significance of the qualitative nature of disorder. This can be done, on the one hand, by 
changing coupling amplitudes with position rather than by cutting some links as performed 
in the present paper, or, on the other hand, by considering topological disorder such as a 
fractal stucture. 
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FIG. 1: The phase diagram for the Miher and Huse map for /i = 3 (top) and /i = 1.9 (bottom). 
The lattice sites are N = L x L with L = 128 and magnetization is computed once T = 2 x 10^ 
iterations have been performed and exponents 7 and /3 are those identified in [l^ . 



FIG. 2: Snapshots for the Miller and Huse CML on the regular lattice for g > Qc- Case j 
g = 0.21 (top) with t = 100 (left), t = 1000 (right). Case /x = 1.9, g = 0.20 (bottom) with t 
(left), t = 1000 (right). The lattice size is L = 128. 
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FIG. 3: The persistence probability for the Miller and Huse maps on the regular lattice for /x = 1.9 
{g = 0.16 is the uppermost curve). The lattice size is L = 128. 




FIG. 4: Phase diagram for Miller and Huse maps for = 3 and dilution: p = 0.05 (top) and 
p = 0.065 (bottom). Lattice size is L = 128. Magnetization is computed once T = 2 x 10^ 
iterations have been performed and exponents 7 and (3 are those identified in 14]. 
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FIG. 5: Critical coupling gc as a function of the dilution parameter p. Miller-Huse model for 
The solid line represents a best fit result yielding ~ O-H for gc = 0.25. 
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FIG. 6: Persistence for the Miller and Huse map for n = 1.9 and dilution p = 0.05. 
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FIG. 7: Snapshots for the Miller and Huse map with dilution (p = 0.05) for g > gdp) in the case 
^ = 3, 5 = 0.23, t = 1000 (top), and in the case n = 1.9, g = 0.20, t = 1000 (bottom). Here the 
size is L = 128. 
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FIG. 8: Various rcscalcd patterns at different p but constant pL^. Parameters are g = 0.2 and 
H = 1.9. Prom top to bottom, L = 128 and p = 0.05; L = 256 and p = 0.0125; L = 512 and 
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FIG. 9: The logistic CML on a regular lattice : the spatial average C{t) = y J2x ^x{t) is shown 
as a function of /j, for the six consecutive times r + l,r + 2,r + 3,r + 4,r + 5,r + 6 where T = 10^ 
for L = 128 and T = 10^ for L = 1024. System sizes and couplings are g = 0.2, L = 128 (top left), 
g = 0.2, L = 1024 (top right), g = 0.1, L = 128 (bottom left), g = 0.1, L = 1024 (bottom right) 
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FIG. 10: (Color online) Snapshots of spin ax{t) (left) and field Cx{t) (right) for fj. = 1.66, g = 0.1, 
L = 64, T = 10^. Two different random initial conditions have been used for rows 1,2 and rows 
3,4 respectively. Row 1 and 3 correspond to time T and rows 2 and 4 to time T + 1. 
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FIG. 11: Probability distribution function (PDF) of variable C{t) over the lattice for the logistic 
CML on a regular lattice (p = 0) with g = 0.1 and fi = 1.66. Computations are performed with a 
size L = 64 with top and bottom figures differing by the random initial condition. Various curves 
refer to six consecutives times T + l, T + 2, T + 3, T + A, T + 5, T + 6 with T = 10^. 
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FIG. 12: Phase diagram for the undiluted logistic map: order parameter ip vs. coupling g (see 
text). Here size is L = 1024, = 1.7, T = 10^. 
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FIG. 13: (Color online) Phase diagram for the logistic CML with coupling g = 0.2, size L = 128 
with binary disorder p = 0.05 and p = 0.1. The ordered case p = is shown as a comparison. 
Spatial average C{t) = y Cx{t) is displayed as a function of fx for six consecutive times T + 
1, T + 2, T + 3, T + 4, T + 5, T + 6 with T = 4.5 x 10^ 
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FIG. 14: The first 100 Lyapunov exponents vs. the normahzed index x = ijl? for the two- 
dimensional Miller-Huse CML with size L = 128. Full (resp. dotted) line refers to a regular lattice 
at ^ = 3 and g = 0.20 < Qc (resp. (7 = 0.21 > gc). The broken line refers to a diluted case g = 0.20 
and p = 0.1. 




FIG. 15: (Color online) The first 100 Lyapunov exponents vs. the normalized index x = i/L'^ for 
the logistic CML. Here the size is L = 128, coupling g = 0.2 and two parameters = 1.65, 1.85 are 
considered. Two curves refer to the the undiluted (p = 0) and diluted (p = 0.1) logistic CML. 



